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We reanalyze the Hubbard-I approximation by showing that it is equivalent to an effective Hamil- 
tonian describing Fermionic charge fluctuations, which can be solved by Bogoliubov transformation. 
As the most important correction in the limit of large U and weak spin correlations we augment this 
Hamiltonian by further effective particles, which describe composite objects of a Fermionic charge 
fluctuation and a spin-, density- or t) excitation. The scheme is valid for positive and negative U . We 
present results for the single particle Green's function for the two-dimensional Hubbard model with 
and without t' and t" terms, and compare to Quantum Monte-Carlo (QMC) results for the param- 
agnetic phase. The overall agreement is significantly improved over the conventional Hubbard-I or 
two-pole approximation. 
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I. INTRODUCTION 

The Hubbard model is the simplest system which pre- 
sumably incorporates the key features of the strong cor- 
relation limit. Understanding this model will be crucial 
for making any progress with cuprate superconductors, 
colossal magneto resistance systems or Heavy Fermions. 
The special problem in this model is that near half- 
filling it represents a dense system of strongly interacting 
Fermions, a situation in which a perturbation expansion 
in U may not be expected to give any meaningful re- 
sults. In dealing with this model, one can then pursue 
two opposing strategies: one might expect that despite 
the strong interaction a perturbation expansion in U re- 
mains a meaningful approximation, and apply conven- 
tional many-body theory. The latter means that one is 
treating the kinetic energy exactly, and results in the va- 
lidity of the Luttinger theorem. On the other hand, by 
adiabatic continuity this approach will never produce an 
insulator at half-filling so that we can be sure of its break- 
down in the limit of large U jt. 
The opposite point of view was taken by Hubbard 
in his approximations, the interaction part of H is treated 
exactly and approximations are made to the kinetic en- 
ergy. This results in the breakdown of the Luttinger theo- 
rem, because the physical electron is effectively split into 
two particles, one of them corresponding to an electron 
moving between empty sites, the other to an electron 
moving between sites occupied by an electron of opposite 
spin. The energies of formation of these effective particles 
differ by [/, whence for large U the single free-electron 
band splits up into two bands formed predominantly by 
the two types of effective particles. Hubbard's approx- 
imations, and related schemes like the so-called 2-pole 
approximations |^-^, have been dismissed by some au- 
thors as unphysical, because they do in fact violate the 
Luttinger theorem away from half-filling. However, in 
a recent QMC study for the paramagnetic phase of the 



Hubbard model |j] we have shown that such criticism is 
entirely unwarranted: the Fermi surface, if measured in 
an 'operational way' from the Fermi energy crossings of 
the quasiparticle band, indeed does violate the Luttinger 
theorem. The doping dependence of the Fermi surface 
volume is qualitatively consistent with Hubbard's results 
and the main discrepancy being the fact that one can 
rather clearly distinguish 4 'bands' in the spectral func- 
tions rather than the 2 bands predicted by the Hubbard- 
I approximation. More generally, exact diagonalization 
studies Q tend to produce for example a photoemission 
spectrum consisting of a relatively narrow 'quasiparti- 
cle band' (which forms the first ionization states) and 
an 'incoherent continuum'. A simple two-band form of 
the spectrum, as produced by the Hubbard-I approxima- 
tion, therefore cannot give a quantitative description of 
the spectrum. Motivated by these numerical results we 
have re-examined the Hubbard I approximation and at- 
tempted to find the most important corrections to this 
scheme for the limit of large U and weak spin correla- 
tions. We will see that the 4-band structure observed in 
the QMC simulations can be reproduced quite well by 
adding two new effective particles which correspond to 
a composite object of a 'Hubbard quasiparticle' and a 
spin-, charge-, or ry-excitation. These composite objects 
actually are the quite obvious leading correction over the 
Hubbard-I approximation, and we will see that includ- 
ing them into the equations of motion leads to an almost 
quantitative agreement with the numerical results for a 
variety of different systems. 



II. A REFORMULATION OF THE HUBBARD I 
APPROXIMATION 



We consider the Hubbard Hamiltonian in particle-hole 
symmetric form, H = Ht + Hjj with: 
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(1) 



Here denotes summation over all pairs of nearest 

neighbors and ^ = c^^ ^Ci^„. 

For bipartite lattices one can make use of two symmetry 
transformations. The first one is the particle-hole trans- 
formation: ^ e^'^'^'c^^ g., where Q = (tt, tt, . . . , tt). If 
the kinetic energy contains only nearest neighbor hop- 
ping, this transformation leaves the Hamiltonian in- 
variant and exchanges the electron addition spectrum 
for momentum k and removal spectrum for momen- 
tum fc -t- Q at half-filling. Similarly the transformation 

| ^ e''5--R«cJ j^, Ci| ^ Cj| inverts the sign of Hu 
At half-filling this allows to transform solutions of the 
positive-U model into those of the negative-U model. 
This transformation implies that the single-particle spec- 
tral function at half-filling is identical for positive and 
negative U . In the appendix some of the operators intro- 
duced in this work and their transformation properties 
under these transformations are listed. 
We proceed with the calculation of the single-particle 
Green's function. As a first step, following Hubbard [Q, 
we split the electron annihilation operator into the two 
eigenoperators of the interaction part: 



These obey [d^^a,Hu] = ^^di^„ and [ci^„,Hu] 



(2) 
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Next, we consider the commutators of these 'effective 
particles' with the kinetic energy. After some algebra, 
thereby using the identity ni_„ = ^^ -\- aSf, we find: 



jeN{i) 



jeN{i} 



(3) 



Here N{i) denotes the z nearest neighbors of site i. Keep- 
ing only the first term in each of the square brackets on 
the r.h.s. (as we will do for the remainder of this section) 
reproduces the Hubbard I approximation. We specialize 
to half-filling {{rii.c) = 1/2) and introduce the Green's 
functions 



(4) 



where a,P € {c, d}. Then, using the anticommuta- 
tor relations {dl^cr^^i^a} = ?^^<T, {cl^cnCi,a} = (1 - riis), 



{d|^,Ci,cr} = {cl^,di^cr} = we obtain the following 
equations of motion: 



idtGc,c = 


^Sit) + ^—Gc^c + 
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(5) 



Taking into account that the ordinary electron Green's 
function G and the Green's function called F by Hubbard 
[|l| can be written as 



G - Gc^c + G- j + Gj^ ~ + G^j, 
r = Gi ^ + 



(6) 



the resulting equations of motions are precisely those de- 
rived in the Hubbard-I approximation: 



idtG^S{t) + {ek-^)G + Ur, 
idtT = ^i6{t)+ekG + Ur). 



(7) 



The present formulation, on the other hand, allows for an 
appealing physical interpretation of the Hubbard-I ap- 
proximation: we introduce free Fermion operators h^, ^ 

and (i^^, which correspond to 'holes' and 'double oc- 
cupancies'. The Hubbard operators are identified with 
these as follows: 



(8) 



Then, we can formally obtain the set of equations of mo- 
tion from the following Hamiltonian for the holes and 
double occupancies: 



2 ^k,a h.a C2 



hk.a^k.a 



(9) 



The Hamiltonian is a quadratic form and readily 
solved by Bogoliubov transformation: 



7- 



t 



7+,fe,tT = -Wfc4^^ + Mfc/iLfc,ff, (10) 

to yield the familiar dispersion relation 



i?±(fe) = i [6fc±y/e| + C/2]. (11) 

To compute the spectral weight of the two bands we use 
Ck,a = -^idk,a + h^-k,s)^ whence 

Z±{k) = ]^{ukTVkf- 

= ifl±_!^V (12) 

Again, this is the correct Hubbard-I result. The above 
discussion shows the physical content of the Hubbard-I 
approximation: the Hamiltonian (^) describes Fermionic 
particle-like and hole-like charge fluctuations, created by 
d\ and h}_i^ g, respectively. These 'live' in a background 
of singly occupied sites. Particle-like and hole- like charge 
fluctuations are created in pairs on nearest neighbors, 
and individually can hop between nearest neighbors. The 
hopping integral for the hole-like particle has opposite 
sign as that for the electron-like fluctuation as it has to 
be, and the hopping integrals for both effective particles 
arc 1/2 times that for the ordinary electrons: this reflects 
the fact that due to the Pauli principle (say) a spin up 
electron added to a 'background' of singly occupied sites 
can propagate to a neighboring site only with a probabil- 
ity of 1/2 (we are neglecting any spin correlations of the 
background of singly occupied sites). The factor of l/-\/2 
in (^) is due to the fact that {cl ^c^ ^) — 1/2. Finally, the 
particle which stands for the double occupancy has an en- 
ergy of formation of L//2, the hole-like particle has an en- 
ergy of — C//2. As already mentioned, the Hubbard I ap- 
proximation therefore describes the splitting of the phys- 
ical electron into two new effective particles which carry 
with them information on the 'environment' in which the 
electron has been created. One of them (dk^a) moves be- 
tween sites occupied by an electron of opposite spin, the 
other one {h^k,a) between empty sites. This is a quite 
appealing physical idea, but the above formulation also 
very clearly highlights the weak points of the Hubbard- 
I approximation. In addition to the mere truncation of 
the commutators (^, which is an uncontrolled approx- 
imation, these are the following: adopting this picture 
we would have to assume that the states hl-^d^ ||0) and 

/ij^^d^ Jo) are distinguishable (and in fact orthogonal to 
one another). This, however, is in general not the case: 
both states have one double occupancy on site j and a 
hole on site i, and the only difference is that they have 
been created in different ways. In fact, using (^ we find 
for their overlap 

^~1{S..-S,), (13) 

where we have assumed a rotationally invariant ground 
state in the second line. The Hubbard-like approximation 



scheme thus should work only for a state with vanishing 
spin correlations - we will therefore henceforth assume 
the spin correlation function {Si ■ Sj) to be zero. 
A second problem is, that the effective Fermions have to 
obey a kind of hard-core constraint - a site cannot be si- 
multaneously occupied by a hole and a double-occupancy. 
This constraint is not accounted for in the derivation of 
the Hubbard-I approximation: the equations of motion 
are obtained from the Hamiltonian (^) by treating the h 
and d as ordinary free Fermion operators. This problem 
is the source of the violation of certain sum-rules when 
the Hubbard-I approximation is applied in the doped 
case, see the discussion given by Avella et al. [|o|. 
One last remark is that the commutators (|^) are invari- 
ant under the particle-hole transformation, i.e. they are 
transformed into each others Hermitian conjugate. This 
remains true for the truncated commutators, which give 
the Hubbard-I approximation, so that the spectral func- 
tion obtained from this is particle-hole symmetric. This 
can also be verified directly using (^l]) and (|2|) . Finally, 
the spectral function is manifestly invariant under sign 
change of /7, as it has to be. 

III. EXTENSION OF THE HUBBARD I 
APPROXIMATION 

We now want to try and derive an improved version 
of the Hubbard-I approximation. Thereby we will ad- 
dress neither the problem of nonorthogonality of different 
hole/double-occupancy configurations nor the hard-core 
constraint - instead, in this work we will restrict ourselves 
entirely to an approximate way of treating the omit- 
ted terms in the commutator relations (I). We expect 
that the present approximation is reasonable for large U 
(where the density of holes/double occupancies is small 
whence the hard-core constraint is of little importance) 
and weak spin correlations (where the nonorthogonality 
problem is small). Throughout we stick to the case of 
half-filling and no spin polarization, {ni^^) = 1/2. We re- 
turn to the basic commutator relations (g|) and consider 
the terms on the r.h.s. which are omitted in the Hubbard- 
I scheme. The second term in the square bracket is 
the Clebsch-Gordan contraction of the spin- 1/2 opera- 
tor Cj^a- and the spin-1 operator Si into yet another spin- 
1/2 operator - it describes the coupling of the created 
hole/annihilated double occupancy to spin excitations. 
This term may be expected to be the most important one 
in the limit of large positive U. The third term describes 
in an analogous way the coupling to density fluctuations, 
whereas the last term describes the coupling to the so- 
called 77-pair excitation [ pT| . One may expect that in the 
case of negative U the last two terms are the important 
ones. 

In keeping with the basic idea of the Hubbard approxi- 
mations, namely to treat the dominant interaction terms 



exactly, we split also the composite operator into eigen- 
operators of the J7-term and define: 
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where hi = rii — (n). Under the positive/negative U 
transformation we have for example Cj,^Sf + Cj,iS^ — > 
■^hidj^^ + Ci,iCi^^Cj ^, i.e. keeping the at first sight unim- 
portant (for positive U) terms involving density and pair- 
ing fluctuations is crucial for maintaining the exact sym- 
metry under sign change of U. We then have 
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(15) 



The operators Cij^a and Di j ^j may be thought of de- 
scribing 'composite objects' consisting of a charge fluc- 
tuation and a spin-, density- or 77-excitation on a nearest 
neighbor. Ultimately these composite operators carry the 
quantum numbers of a single electron, i.e. spin 1/2 and 
charge 1. We also note that under particle-hole transfor- 
mation 



C, 



(16) 



i.e. the composite particles transform in the same way 
as the Ci^CT and di^. Moreover, the commutators ( [l5| ) 
are invariant under particle-hole transformation, i.e. this 
transformation transforms them into each others Hermi- 
tian conjugates. 

We now enlarge the set of Green's functions by allow- 
ing a,/? e {c,d,C,D} in (^; in the language of the 
'effective Fermions' this means that we are introducing 
additional Fermions corresponding to the composite ob- 
jects. To obtain a closed system of equations of motion 
for these Green's functions, we need equations for the 
and Gj^-. As a first step we turn to the commutators 
[Ci,j,a,Ht\ and [Dij^a, Ht]. Here we have to distinguish 
three cases (see Figure |l|) : 

a) the hopping term may act along the bond and 
transport the hole back from j to i. 

b) it may transport the hole even further away from i 

c) it may transport the spin- density- or 77-excitation 
away from site i. 

If we want to restrict ourselves to the 4 types of opera- 
tors Cij,a-, Di j^a, Ci,a- and di^a-, we have to neglect the 
contributions from the processes b and c. These pro- 
cesses would produce 'strings' of excitations along a path 



of length 2 lattice spacings, and we would have to in- 
troduce even more complicated operator products to de- 
scribe these. Restricting ourselves to processes of the 
the type a, i.e. replacing Ht Z]o-(4,o-Cj,o- + H.c), 

straightforward computation gives the surprisingly 
simple result 

3t ~ t - Zt 

it - t " Zt 

[AjT>^^*] = 2"^^'*'T + 2^^'^'^ " T^'^^'^ " ^^^^ 

Again, these relations are particle-hole invariant, i.e. 
they are turned into each others Hermitian conjugates 
by particle-hole transformation. In passing we note that 
had we reduced the operators A.jT a^*i ^jjT ^'^ comprise 
only the terms involving spin excitations (as might seem 
appropriate in the case of large positive U) , the commu- 
tators would have been much more complicated and in 
fact the 'Hamilton matrix' Hk to be defined below would 
have been non-Hermitian. 
Next, we need the anticommutators 



(18) 



Taking the expectation value in the ground state most of 
the terms vanish on the basis of symmetries: Cj^j^d^ -I- 

vanishes due to inversion symmetry of the ground 

state, Cj_|c| I — and hi vanish due to particle-hole 

symmetry at half-filling. All terms containing unpaired 
spin operators vanish if we assume that the ground state 
is invariant under spin rotations (which excludes ferro- 
or antiferromagnetic solutions). Finally we obtain: 



({A,„t,£[t}) = 



(19) 



It is is easy to see that the expressions whose expectation 
values are taken are invariant under particle-hole trans- 
formation, and under the positive/negative U transfor- 
mation we have: 



rLiTij 



(20) 



The expectation value of the anticommutator would be 
invariant under this positive/negative-U symmetry. 



For large positive U the terms and ] |C] ^Cj^^Cj^i have 
a neghgible expectation value and the only important 
term comes from the spin correlation. In keeping with 
our above remarks concerning the role of spin correlations 
in the Hubbard I approximation we will henceforth take 
the r.h.s. of ( p^ to be zero. As was discussed above, the 
Hubbard-I approximation implicitly assumes (Si-Sj) = 0, 
and we will therefore keep this value also in . We will 
discuss the consequences of not making this approxima- 
tion later on. 

Using the above commutators and (expectation values 
of) anticommutators we are now in a position to set up a 
closed system of equations of motion. In the following we 
give explicit formulas only for a ID chain with nearest- 
neighbor hopping, but the generalization to higher di- 
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mensions and/or longer range hopping integrals will be 
self-evident. We introduce the Fourier transforms 




3 



(and analogously for the i)'s) and define the vector 
Gc = (Ccc, Gjg,G'(j,^ g,G^_ g,G^^ g,G'£,__g). Here C± 

is shorthand for C±^a{k). Combining (||), (|lj) and 
( |l7| ) , and performing a spatial Fourier transformation the 
equations of motion are readily found to be 

{idt - Hk)Gc - S{t)B, (21) 

where the Hermitian matrix Hk is given by 



_^gifcx/2 




^ie'^j2 


\ 




^g — ifcx/2 


fgifex/2 


t 





3t 




2 ' 
U 


3t 


2 






2 ' 


2 ' 




3t 


U 


t 




2 ' 


2 ' 


2 







t 

2 ' 


U 
2 


/ 



_ig-«fe./2 

^g-ifcx/2 

_U_ 
2 

t_ 

2 



2 



and Be = O7 0, 0, 0). The equation system (^ ) can 
be solved for each momentum and frequency by using 
the spectral resolution of the Hamilton matrix Hk- This 
yields the Green's functions Ggg and Gjg for each mo- 
mentum and frequency. In an analogous way we can also 
derive an equation system for G^j and Gjj. Thereby 
the matrix Hk stays unchanged, whereas the r.h.s. is 
changed into Bd = (0,1/2,0,0,0,0). Finally, the full 
electron Green's function is obtained by adding up the 
four Green's functions according to (^). Upon forming 
the Laplace transform G(fe, z) we finally obtain the spec- 
tral density 



A{k, u}) = - lim -'^G{k, uj + iij). 

77^0 TT 



(22) 



In passing we note that this way of computing A{k,uj) 
guarantees the validity of the sum rule 



A{k,u!)dijj = 1. 



(23) 



Since the particle-hole symmetry of the relations (15), 
( p^ ) and (|l^) in turn guarantees particle-hole symmetry 
of the entire spectral function, we find that the sum-rule 
for the particle number is fulfilled automatically: 



E 



A{k, uj)dtJ 
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(24) 



As a last remark we note that going over to higher di- 
mensions or adding longer ranged hopping integrals poses 
no problem whatsoever - for each spatial dimension a we 



have to add four additional rows and columns containing 
the Cij and Di^ where i and j are nearest neighbors 
in ±a-direction. Similarly, if we add an additional hop- 
ping integral t' between 2""^ or 3'''^ nearest neighbors (the 
number of whom we denote by z') to the Hamiltonian, 
we have to add 2z' rows and columns, containing the Cij 
and Di j with 2'"* or 3'''' nearest neighbors i and j. In 
each case these additional rows and columns contain only 
mixing terms amongst themselves or with Gc,c and G^~, 
so that the extension is completely trivial. 



IV. COMPARISON WITH NUMERICS 

Following the discussion in the preceding section we 
can calculate the full electron Green's function, includ- 
ing the (presumably) most important corrections over the 
Hubbard I approximation in the limit of weak spin cor- 
relations and large U. We now proceed to a comparison 
of the obtained results results for the spectral density 
AlkjUj), with the spectral density obtained from Quan- 
tum Monte-Carlo (QMC) simulations. Thereby the tem- 
perature for the QMC simulation, T = 0.33i, was cho- 
sen such that the spin correlation length is only approxi- 
mately 1.5 lattice spacings - the results thus are probably 
quite representative for the paramagnetic regime which 
our approximation aims to describe. Moreover, the value 
oi U/t — 8 is already rather large, so that we may also 
hope to have a small density of holes/double occupancies 
and the neglect of the hard-core constraint be justified. 
As a general remark concerning the QMC spectra we note 



that the MaxEnt procedure used for the analytic contin- 
uation to the real axis es most reliable for 'features' with 
large weight - this means that the position of tiny peaks 
is less accurate than that of large ones. 
Figure ^ then compares the spectral density obtained 
from the Hubbard I approximation, our extended Hub- 
bard approximation (EHA) and QMC simulation. The 
Hubbard I approximation gives only a relatively crude 
fit to the actual spectral density obtained by QMC. The 
extended Hubbard approximation, on the other hand, 
gives an all in all quite correct description of the spectral 
density. Out of the 10 = 2 + 4 -I- 4 bands produced by 
diagonalizing i/^, only four bands do have an apprecia- 
ble spectral weight. These four intense bands correspond 
rather well to 4 broad 'bands' of intense spectral weight 
which can be roughly identified in the QMC result. It is 
interesting to note that a recent strong-coupling expan- 
sion for the Hubbard model by Pairault et al. |]l3| also 
produced a 4-band structure, although only for the ID 
model. The dispersion of the spectral weight along the 
bands is also reproduced quite satisfactorily be the ex- 
tended Hubbard approximation. The main difference is 
the apparently strongly fc-dependent width of the spectra 
produced by QMC, which however is outside the scope 
of the present approximation, which (in the limit 77 ^ 0) 
produces sharp J-peaks without any broadening. A more 
severe deficiency of the EHA is, that it tends to predict 
too high excitation energies, resulting in a somewhat too 
large value of the Hubbard gap. In any way, however, the 
magnitude of the Hubbard gap comes out better than in 
the Hubbard-I approximation. 

We proceed to the Hubbard model with an additional 
hopping integral t' between 2'"^ nearest (i.e. (l,l)-like) 
neighbors. Figure ^ again compares the Hubbard I ap- 
proximation, the EHA, and the results of a QMC simu- 
lation on an 8 X 8 lattice. The agreement between EHA 
and the QMC result is again quite good, the main dis- 
crepancy being again an overall overestimation of the 
binding energies. On the other hand, the apparent 4- 
band structure, the dispersion of the peak energies and 
the spectral weight agrees well with the numerical result. 
In particular, the 2 bands in inverse photoemission (i.e. 
a; > 0) predicted by the EHA can be seen very clearly 
in the QMC spectra). All in all the agreement is even 
better than for the case t' = 0, which most probably is 
due to the fact that the spin correlations are weaker with 
i', so that the assumption of an uncorrelated spin state 
is better justified in this case. 

We proceed to the case of a hopping integral integral t" 
between 3'"'' nearest (i.e. (2,0)-like) neighbors. Here we 
have chosen t" /t = 0.25, because for the larger value 
t" /t = 0.5 the QMC simulation still predicted a metallic 
state at U/t = 8. Figure ^ shows the spectral functions 
for t" /t — 0.25. Again, we can roughly identify 4 bands 
and there is reasonable agreement for the dispersion. The 
weight of the quasiparticle band in photoemission near 



(tt, tt) (at uj « —2) is not reproduced very well by the 
EHA, but again the ubiquitous 4-band structure is rather 
clearly visible. 

Next, we turn to a somewhat indirect check of the ap- 
proximation. The ordinary electron creation operator is 
the symmetric combination of the Hubbard operators. 
However, we might also define the antisymmetric combi- 
nation 



.(1 



(25) 



This operator has the same quantum numbers as the orig- 
inal electron operator, and therefore obeys the same se- 
lection rules. It follows that when acting on the ground 
state, this operator probes the same final states as the 
electron operator, the only difference being the matrix 
element viz. the spectral weight of the respective peak 
in the spectral density. In fact the Green's function 



can be expressed as 



G — Gcc ~ Gdc ~ G, 
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(27) 



It therefore is easy to calculate within our approximation, 
and comparison with the QMC result provides an addi- 
tional check for our description of the electronic struc- 
ture. Note that the operator Ck.a enhances the practi- 
cally dispersionless bands at w = ±6 - since these bands 
now have a larger weight, their position and dispersion 
are more reliable than in the 'ordinary' photoemission 
spectrum. This is shown in Figure |^, where indeed quite 
good agreement is found between the EHA and the nu- 
merical result. 

Finally, we turn to the discussion of choosing a nonvan- 
ishing expectation value of the anticommutator in (|is|), 
i.e. we assume that ({Ci j,cr, c|^}) = Sjjx 7^ 0. This 
changes the rhs of the equation system ( pi|) to Be = 

(1/2, 0, \J^Xj . . .) but leaves the matrix Hk unchanged. 
In other words, the dispersion of the bands, which is de- 
termined by the eigenvalues of stays unchanged and 
only the spectral weight of the peaks changes. Moreover 
the sum rules (23) and (^J) retain their validity also in 
this case. For large positive U and n = \ charge fluctua- 
tions will be strongly suppressed and double occupancies 
will have a small probability, so that the dominant con- 
tribution to X comes from the spin correlation function 
{Si ■ Sj), whence we should choose x < 0. Assuming for 
example a negative x of moderate value, x — —0.2, then 
leads to little change in the calculated spectral density 
(see Figure ^): the same two bands which had a large 
spectral weight for a; = retain a large spectral weight 
also in this case. There is, however, a rather undesir- 
able feature associated with the bands with small spec- 
tral weight: numerical evaluation shows, that for some 



regions of fc-space these bands acquire a small but nega- 
tive weight. 

The physical origin of this problem is probably related 
to nonorthogonalities of basis states: in principle we 
could interpret the matrix Hk as a Hamiltonian describ- 
ing (in 2D with only nearest neighbor hopping) 6 types 
of Fermionic 'effective particles'. Quite generally, the 
anticommutator-relation {a, 5^} — x ^ implies that 
the wave functions corresponding to the Fermi operators 
and 6^ are non-orthogonal. While an exact overlap 
matrix can never have negative eigenvalues but at most 
develop zero eigenvalues (indicating that the set of basis 
states is overcomplete), any approximation to the matrix 
elements may lead, as an artefact of the approximation, 
to negative eigenvalues. Since we are using only approx- 
imate values for the overlaps, it may happen that we ob- 
tain states with a nominally negative norm, whence we 
can get poles of negative weight. Setting x = through- 
out amounts to assuming that all our effective particles 
are orthogonal to one another and obviously removes the 
problem with nonorthogonalities. This seems reasonable, 
because we are neglecting overlap terms proportional to 
the spin correlation function already in the Hubbard-I 
approximation. The lesson then is basically the same as 
discussed already before: the Hubbard-I approximation 
is well defined only when applied to an 'ideally param- 
agnetic' state with no correlations of finite range, and 
applying it to a state with finite spin correlations repre- 
sents an approximation. 

V. CONCLUSION 

In summary, we have investigated the most impor- 
tant corrections over the Hubbard-I approximation in 
the limit U /t oo and electron density n = 1. We 
have seen that the Hubbard-I approximation describes 
charge fluctuations on a 'background' of singly occupied 
sites, which is moreover assumed to have zero spin corre- 
lations. The charge fluctuations are point-like, and corre- 
spond to an electron moving between empty sites and an 
electron moving between singly occupied sites. We note 
that a very similar construction can also be applied to 
the Kondo lattice and in fact reproduces the single 
particle spectra very well. This is probably due to the 
fact that the Kondo lattice has a unique ground state in 
the limit of zero kinetic energy, whereas the ground state 
of the Hubbard model is highly degenerate in the case 
t = 0. 

In our extended scheme for the Hubbard model we have 
augmented these point-like charge fluctuations by addi- 
tional 'particles' which arc composite in character and 
consist of a point-like charge fluctuation together with a 
spin-, density-, or 77-excitation. Comparison of the ob- 
tained single-particle spectral density with QMC results 
for a variety of systems showed a quite reasonable agree- 



ment. In particular the apparent 4-band structure seen 
in the numerical spectra finds its natural explanation in 
the extended Hubbard approximation. We also note that 
QMC simulations where the spectra of the composite ex- 
citations have actually been computed , further sup- 
port the present interpretation. We thus have a quite suc- 
cessful method of computing the full quasiparticle band 
structure of the Hubbard model, at least in the param- 
agnetic case and at half-filling. 

The present scenario for the nature of the composite exci- 
tations also allows to make a connection with various the- 
ories for the hole motion in an antiferromagnet [p^|-pl|. 
There, one is describing holes dressed by antiferromag- 
netic spin fluctuations. When acting on the Neel state 
the operators Cij^cr obviously describe precisely a hole 
together with a 'spin wave' on a nearest neighbor or, put 
another way, a 'string' of length one. The terms which 
were omitted from the equation of motion for the Cij^a 
then would correspond to strings of length two and so 
on. While such longer-ranged strings are apparently of 
minor importance in the paramagnetic phase, one may 
expect that they become more and more important for 
the description of the dispersion the stronger the antifer- 
romagnetic correlations. The relative importance of such 
longer ranged strings therefore may be the mechanism for 
the crossover from the Hubbard-I like dispersion in the 
paramagnetic phase to the spin-density-wave like disper- 
sion in the antiferromagnetic phase. Similarly, one might 
think of formulating the entire Hubbard-I approximation 
also in the antiferromagnetic phase, by constructing the 
Hamiltonian for charge fluctuations explicitly for a Neel 
ordered spin background. 
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VI. APPENDIX 

Transformation properties of various operators under 
particle- hole and positive/negative U transformation. 
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FIG. 2. Single particle spectral function for the Hubbard 
Model with U/t = 8 from the Hubbard I approximation, the 
extended Hubbard approximation, and QMC simulations on 
a 20 X 20 lattice at temperature T = t/3. In this as well as in 
all following figure, the approximate spectra have been given 
an artificial Lorentzian broadening ri = 0.20t. To compen- 
sate for the stronger broadening the QMC spectra have been 
multiplied by an additional factor of 2. 
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FIG. 1. Possible hopping processes which couple to Ci.j,a- 
The cross denotes the spin-, density- or 7;-excitation, the black 
dot the hole. 
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FIG. 3. Single particle spectral function for the Hubbard 
Model with U/t = 8, t' = t/2 from the Hubbard I approxima- 
tion, the extended Hubbard approximation, and QMC simu- 
lations on an 8 X 8 lattice at temperature T = t/3. To com- 
pensate for the stronger broadening the QMC spectra have 
been multiplied by an additional factor of 4. 
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FIG. 4. Single particle spectral function for the Hubbaxd 
Model with U/t = 8,t" = t/4 from the Hubbard I approxima- 
tion, the extended Hubbard approximation, and QMC simu- 
lations on an 8 X 8 lattice at temperature T = t/3. To com- 
pensate for the stronger broadening the QMC spectra have 
been multiplied by an additional factor of 4. 
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FIG. 5. Spectral function of the c-operator for the Hub- 
bard Model with U/t = 8 from the Hubbard I approximation, 

the extended Hubbard approximation, and QMC simulations 
on an 8 X 8 lattice at temperature T = t/3. To compen- 
sate for the stronger broadening the QMC spectra have been 
multiplied by an additional factor of 4. 
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FIG. 6. Spectral function calculated with the extended 
Hubbard approximation for x = —0.2. 



